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00 ■ Abstract 

o 

I The purpose of this note is to report predictions of the binding energy of 

the A hyperon in nuclear matter using the latest version of the Jiilich nucleon- 
C . hyperon meson-exchange potential. Results from a conventional Brueckner cal- 

culation are compared with previously reported values. A calculation including 
\^ ■ Dirac effects on the A single-particle potential is also presented. Issues encoun- 

tered in Dirac calculations with nucleon-hyperon potentials are discussed. 
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There are important motivations for including strange baryons in nuclear matter. 



Q\ \ The presence of hyperons in stellar matter tends to soften the equation of state 

^ ' (EoS), with the consequence that the predicted neutron star maximum masses be- 

Q ■ come considerably smaller. With recent constraints allowing maximum masses larger 

than previously accepted limits [1], accurate microscopic calculations which include 
O I strangeness (in addition to other important effects, such as those originating from rel- 

ativity), become especially important and timely. Furthermore, as far as terrestrial 
nuclear physics is concerned, studies of hyperon energies in nuclear matter naturally 
complement our knowledge of hypernuclei. 

Microscopic calculations of nuclear matter properties have been reported ear- 
^ ! lier within the non-relativistic Brueckner-Hartree-Fock framework (see, for instance, 

Refs. [2], [3]), using the Nijmegen [1] and/or the Jiilich [5], [6] nucleon-hyperon (NY) 
meson-exchange potentials. 

In this work we use the most recent meson-exchange NY potential from the Jiilich 
group [7j. Given that there are significant differences between this and the previous 
(energy independent) version of the Jiilich NY potential |6|, it will be interesting to 
see how those differences reflect onto G-matrix calculations. Before moving on to our 
final objective, which is a microscopic determination of the EoS for hyperonic matter, 
we will first confront a much simpler scenario, namely the one of nuclear matter 
at some Fermi momentum kp in the presence of a "A impurity". This calculation, 
the outcome of which is the binding energy of a A hyperon in nuclear matter, will 
allow us to explore the behavior of the new NY potential in nuclear matter before 
addressing more involved situations. Some empirical information is available for the 
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A binding energy from analyses of (tt, K) and (iT, tt) reactions and studies of energies 
of hj^ernuclei [8] . 

We will also report on our prediction for the A binding energy from a Dirac- 
Brueckner-Hartree-Fock (DBHF) calculation. To the best of our knowledge, such 
calculation has not been performed before. 

2 Calculations using the non-relativistic Brueck- 
ner G-matrix 

For matter with non-vanishing hyperonic density, the nucleon. A, and S single- 
particle potentials are the solution of a coupled self-consistency problem, which reads, 
schematically 



In the equations above, G^n^ Gny, and Gyy, (X,y = A, S), are the nucleon- 
nucleon, nucleon-hyperon, and hyperon-hyperon G-matrices at some nucleon and hy- 
peron densities defined by the Fermi momenta kp and k^. 

Following the earlier calculation from Ref. [6], which we want to compare with, 
we make the following assumptions and approximations: 

• We consider the case of symmetric nuclear matter at some Fermi momentum 
kp in the presence of a "A impurity", i.e. kp ^ 0. 

• We take the single-nucleon potential from a separate calculation of symmetric 
matter [9]. Notice that the A potential is quite insensitive to the choice of Un, 
as reported in Ref. [6] and as we have observed as well. 

• Both A and S are included in the coupled- channel calculation of the NY G- 
matrix, but free-space energies are used for the latter. 

The parameters of the A potential, on the other hand, are calculated self-consistently 
with the Gny interaction, which is the solution of the Bethe-Goldstone equation with 
one-boson exchange NY potentials. In the Brueckner calculation, density-dependent 
effects come in through angle-averaged Pauli blocking and dispersion. 

For some NY initial state, the starting energy that enters in the scattering equation 

is 




(1) 
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For nucleons and A's, we write the single-particle energy, in the non-relativistic case, 
as ^ 

^'^P^ = + ^^^Pi") + (3) 

{i = N,A). (Angle-averaged) Pauli blocking is applied to all intermediate two-baryon 
states. The integration in the Bethe-Goldstone equation is handled using standard 
methods to eliminate the angular dependence [0]. We adopt the continuous choice 
for the single-particle potential. Once the latter is obtained, the value of —U\{p\) at 
Pa=0 provides the A binding energy in nuclear matter, B\. 

In this work, we will apply the latest meson-exchange model by the Jiilich group 
[7], which will be denoted by NY05. As shown and discussed extensively in Ref. [7], 
there are several remarkable differences between this model and the previous NY 
Jiilich potential [6] (NY94). The main new feature of NY05 is a microscopic model 
of correlated tttt and KK exchange to constrain both the a and p contributions 
[7]. The usual one-boson-exchange contributions from the lowest pseudoscalar and 
vector meson multiplets are also present, with the coupling constants determined 
by SU(6). This makes the long and intermediate range parts of the potential well 
determined. New short-range features are included through the ao(980) meson and 
a strange scalar meson with a mass of approximately 1000 MeV. Both of these are 
parametrized phenomenologically in terms of one-boson exchanges in the respective 
spin-isospin channels. The model describes well the available data on integrated as 
well as differential cross sections [7|. Also, the hypertriton binding energy is well 
reproduced [lOj. 

There are, though, some major quantitative differences between NY05 and NY94 
in specific partial waves. Most noticeably, the new model predicts a considerably 
larger scattering length in the singlet channel. These differences turn out to have a 
large impact on in-medium predictions. With the new model, we obtain considerably 
more attraction than Reuber et al. [B], approximately 50 MeV at A;^=1.35 fm~^ for 
rather than 30 MeV. A value of 49.7 MeV has been reported by H. Polinder with 
NY05 [11]. Notice that a value of 30 MeV is generally accepted as the "empirical" 
one [8], which opens some interesting questions: 

• The value predicted by NY05 for the hypertriton binding energy is 2.27 MeV, 
in satisfactory agreement with the experimental value of 2.354(50) MeV. How 
reliable is the "empirical" value for the A binding energy in nuclear matter? 

• How will the additional attraction impact the EoS and neutron star predictions? 
How will those predictions compare to the most recent constraints? 

To better highlight the potential model dependence of the predicted A binding 
energy, we show in Table 1 how selected partial waves contribute to it. For each 
partial wave, column NY94 shows the results given in Ref. [llj, whereas column 
NY05 are the predictions from this work. Clearly the largest contribution to the 
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Table 1: Contributions to the A binding energy from selected partial waves as obtained 
in non-relativistic Brueckner-Hartee-Fock calculations. NY05 and NY94 indicate the 
Jiilich NY potentials from Ref. ff] and Ref. [6], respectively. 



Partial wave 


NY94, BHF 


NY05, BHF 




3.6 


8.73 




27.2 


37.69 




-0.6 


0.69 




-2.0 


0.13 




0.8 


3.27 


Total (all states) 


29.8 


51.27 



model dependence originates from the iS-waves, although the contribution from the 
P-waves is also dramatically different between the two sets of predictions (but much 
smaller than the one from the S*- waves). 

The NY05 entries agree well with those from Ref. [11] for the same potential; the 
relatively minor differences are possibly due to different choices for the single-particle 
spectrum, the handling of the angular dependence in the scattering equation, and, to 
a very small extent, the choice of the nucleon potential. 

3 Dirac effects on the A binding energy 

The relation between the non-relativistic Brueckner approach and the relativistic 
framework (known as Dirac-Brueckner-Hartree-Fock, DBHF) has been discussed for 
a long time. Already in Ref. [12] it was shown how relativistic effects tie in with virtual 
excitations of pair terms. Lately, these concepts have been revisited in more detail 
|13j and with similar conclusions. In short, the Dirac effect on the EoS of nucleonic 
matter is an essential saturating, and strongly density dependent, mechanism, which 
effectively accounts for the class of three-body forces originating from virtual nucleon- 
antinucleon excitations. When hyperon degrees of freedom are included, for reasons 
of consistency, those should then be subjected to the same correction. 

We have incorporated DBHF effects in the present calculation, which amounts to 
involving the A single-particle Dirac wave function in the self-consistent calculation 
through the A effective mass, m^. Similarly to what is done for nucleons |T4j, we fit 
the single-particle energy for A's using the ansatz 

e^{p) = ^{mlY + p^ + U^ (4) 

with m\ = m\ + Ug, and Ug and Uy the scalar and vector potentials of the A baryon. 

A problem with the Jiilich NY potential in conjunction with DBHF calculations 
is the use of the pseudoscalar coupling for the interactions of pseudoscalar mesons 
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(pions and kaons) with nucleons and hyperons. For the reasons mentioned above 
(that is, the close relationship between Dirac effects and "Z-diagram" contributions), 
this relativistic correction is known to become unreasonably large when applied to 
a vertex involving pseudoscalar coupling. On the other hand, the gradient (pseu- 
dovector) coupling (also supported by chiral symmetry arguments) largely suppresses 
antiparticle contributions. To resolve this problem, one can make use of the on-shell 
equivalence between the pseudoscalar and the pseudovector couphng, which amounts 
to relating the coupling constants as follows: 

_ rrij + rrij 

9ps Jpv ) 

where gps denotes the pseudoscalar coupling constant and /p„ the pseudovector one; 
mpsi rrii, and mj are the masses of the pseudoscalar meson and the two baryons in- 
volved in the vertex. This procedure can be made plausible by writing down the 
appropriate one-boson-exchange amplitudes and observing that, redefining the cou- 
pling constants as above, we have (see Ref. [H] for the two-nucleon case) 

Vp, = Vp, + (6) 

where the ellipsis stands for off-shell contributions. Thus, the pseudoscalar coupling 
can be interpreted as pseudovector coupling where the off-shell terms are ignored. 
This is what we apply in our DBHF calculations. 

In the coupled channel calculation, evaluation of Gna involves the transition po- 
tentials Vna^na, Vna^n's, and V^s^ny; (all with total channel isospin equal to 1/2), 
plus the corresponding exchange diagrams. Because in the present scenario (of a A 
impurity in nucleonic matter) the S hyperon is not given an effective mass, Dirac 
effects are applied only in Vna^na- A diagram where not all of the baryon lines are 
Dirac-modified may yield a Dirac effect that is artificially skewed. Moreover, since we 
find that the net contribution to the A binding energy from the coupling to the NT, 
channels is rather small (~ 1.3 MeV), we anticipate Dirac effects from those channels 
to be negligibly small. 

Finally, a comment is in place concerning meson propagators. In standard DBHF 
calculations [H], the so-called Thompson equation (a relativistic three-dimensional 
reduction of the Bethe-Salpeter equation) is used for two-baryon scattering. In the 
Thompson formalism, static propagators are employed for meson exchange, i.e. 

I (7) 

where rria denotes the mass of the exchanged meson and (f are the baryon momenta 
in their center-of-mass frame before and after scattering. The Jiilich NY potentials 
are based upon time-ordered perturbation theory |5j and use a meson propagator 
given by 

—i T 

Ua{z - Ei - Ej - Ua) 
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Table 2: Contributions to the A binding energy from selected partial waves obtained 
with the Jiilich NY05 potential in a DBHF calculation. 



Partial wave 


NY05, DBHF 




8.14 




36.45 




0.07 




-1.19 


'-'P2 + ''F2 


3.21 


Total (all states) 


47.4 



with Ua = + io' ~ Q}"^'-, Ei = + (i"^ and Ej = ^m?^ + <p are baryon energies 

and z is the starting energy of the two-baryon system. In order to eliminate the 
energy dependence, Reuber et al. |6] replaced the original z with 

z = -{nil + m2 + + rrii) , (9) 

where the mj's denote the baryon masses of the four legs in the one- meson exchange 
diagram. In any case, the Jiilich meson propagator involves the baryon masses. Re- 
placement of these free-space masses with in-medium values would create medium 
effects on meson propagation which we do not wish to include in our nuclear matter 
calculations. The reason for keeping free-space masses in the meson propagator is 
twofold. First, standard DBHF calculations do not include medium effects on me- 
son propagation as they typically use Eq. (7), which does not not depend on baryon 
masses. Second, medium effects on meson propagation constitute a separate class of 
effects that we are not concerned with in the present context and are typically not 
perceived as part of the DBHF approach. 

Having taken the steps described above, we proceed to the DBHF calculation and 
find a moderate reduction of i?A, by approximately 4 MeV, due to Dirac effects. This 
is roughly 50% of the corresponding effect on the nucleon potential. In Table II, 
contributions from selected partial waves are again shown. A large part of the effect 
can be attributed to increased repulsion in S and P waves, especially ^Pi. 

4 Conclusions 

We have reported on non-relativistic and Dirac-Brueckner-Hartree-Fock predictions 
of the A binding energy in nuclear matter at normal nuclear density. The magnitude 
of the Dirac effect is approximately 1 /2 of the corresponding effect on the binding of 
a nucleon in nuclear matter. 
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We also noticed and discussed the remarkably different predictions of this "ob- 
servable" as obtained using the 2005 or the 1994 versions of the Jiilich NY potentials. 
Given that there are noticeable differences between the two free-space potentials, 
some potential model dependence is to be expected. In this case, though, it would 
be appropriate to say that those free-space differences are considerably "amplified" 
in the nuclear matter calculation. Notice that this comparison was done between 
predictions as obtained from conventional (BHF) calculations. But the conclusions 
would be unimpacted by the Dirac effect, which is much smaller than the differences 
originating from the use of the two potential models. 

The natural extension of this preliminary calculation will be a fully self-consistent 
DBHF calculation of Un, Ua, and t/s for diverse N, A and S concentrations. Such 
work is in progress. 
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